Strong-coupling approach for strongly correlated electron systems 
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A perturbation theory scheme in terms of electron hopping, which is based on the Wick theorem 
for Hubbard operators, is developed. Diagrammatic series contain single-site vertices connected 
by hopping lines and it is shown that for each vertex the problem splits into the subspaces with 
"vacuum states" determined by the diagonal Hubbard operators and only excitations around these 
vacuum states are allowed. The rules to construct diagrams are proposed. In the limit of infinite 
spatial dimensions the total auxiliary single-site problem exactly splits into subspaces that allows 
to build an analytical thermodynamically consistent approach for a Hubbard model. Some ana- 
lytical results are given for the simple approximations when the two-pole (alloy-analogy solution) 
and four-pole (Hartree-Fock approximation) structure for Green's function is obtained. Two poles 
describe contribution from the Fermi-liquid component, which is dominant for small electron and 
hole concentrations ("overdoped case" of high-Tc's), whereas other two describe contribution from 
the non- Fermi liquid and are dominant close to half-filling ("underdoped case"). 
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I. INTRODUCTION 

Many unconventional properties (e.g., metal-insulator 
transition, electronic (anti)ferromagnetism) of the nar- 
row-band systems (transition metals and their com- 
pounds, some organic systems, high- Tc superconductors, 
etc.) can be explained only by the proper treatment of 
the strong local electron correlations. The simplest mod- 
els allowing for the electron correlations are a single-band 
Hubbard model with on-site repulsion U and hopping 
energy t and its strong-coupling limit {U ^ t): t — J 
model. Recent studies of the Hubbard-type models con- 
nected mainly with the theory of high- Tc superconduc- 
tivity and performed in the weak- [U < 4t) and strong- 
{U » t) coupling lirpits, elucidate some important fea- 
tures of these models.El But still a lot of problems remains, 
especially for the U ^ t case where there are no rigorous 
approaches. 

Such approaches can be built using systematic Dertur- 
bation expansion in terms of the electron hoppina -using 
diagrammatic technique for Hubbard operators .u'cl One 
of them was proposed for the Hubbard (U = oo limit) and 
t — J models.Elfl The lack of such approach is connected 
with the concept of a "hierarchy" system for Hubbard 
operators when the form of the diagrammatic series and 
final results strongly depend on the system of the pairing 
priority for Hubbard operators. On the other hand it is 
difficult to generalize it on the case of the arbitrary U. 

In the last decade the essential achievements in the the- 
ory of the strongly correlated electron systems are con- 
nected with the development of the dynamical mean-field 
theory (DMFT) proposed by Metzner and VoUhardtlD 
for the Hubbard model (see also Ref. || and references 
therein). DMFT is a nonperturbative scheme that al- 
lows to project the Hubbard model on the single impu- 
rity Anderson model and is exact in the limit of infinite 
space dimensions (d ^ oo). There arc no restrictions on 



the U value within this theory and it turns out to be 
useful for intermediate coupling {U ~ t) for which it en- 
sures the correct description of the metal-insulator phase 
transition and determines the region of the Fermi-liquid 
behavior of the electron subsystem. Moreover, some 
class of the binary-alloy-type models (e.g., the Falicov- 
KimbalLmodel) can be studied almost analytically within 
DMFT.Ij But in the case of the Hubbard model, the treat- 
ment of the effective single impurity Anderson model is 
very complicated and mainly computer simulations [ex- 
act diagonalization of the finite-sized systems or quan- 
tum Monte Carlo (QMC)] are used, which calls for the 
development of the analytical approaches. 

The first analytical approximation proposed for tka 
Hubbard model was a simple Hubbard-I approximationll3 
(see Ref. ^ for its possible improvement) which is cor- 
rect in the atomic (t = 0) and band {U = 0) lim- 
its but is inconsistent in the intermediate cases and 
cannot describe metaidnsulator transition. Hubbard's 
alloy-analogy solutionllj (so-called Hubbard-III approx- 
imation) incorporates into the theory an electron scat- 
tering on the charge and spin fluctuations that al- 
lows us to give qualitative description of the changes 
of the density-of-state at the metal-insulator transition 
point. Hubbard-I and Hubbard-III approximations in- 
troduces two types of particles (electrons moving be- 
tween empty sites and electrons moving between sites 
occupied by electrons of opposite spin) with the dif- 
ferent energies that differ by U and form two Hub- 
bard bands. Related schemes of the so-called two-pole 
approximationsjljo which are justified by the t/U ^ 1 
perturbation theory expansions axe-jalso considered. 
However, in the recent QMC studieslijO it is clearly dis- 
tinguished four bands in the spectral functions rather 
than the two bands predicted by the two-pole approxi- 
mations. Such four-band structure is reproduced bj^-jthe 
strong-coupling expansion for the Hubbard modeltZI in 
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the one-dimensional case. Within other pajj>«roaches let 
us mention nonpecossing approximationJljO Edwapds-, 
Hertz approach,E3Eil iterative peEtuxbation theoryj23'E2l 
alloy-analoffiji.pbased approaches ,Ot3 and linked cluster 
expansions ,E3EZI which are reliable in certain limits and 
the construction of the pthermodynamically consistent 
theory still remains openBS 

The aim of this paper is to develop for Hubbard-type 
models a rigorous perturbation theory scheme in terms 
of electron hopping that is based on the Wick theorem 
for Hubbard operatora3a and is valid for arbitrary value 
oi U {U < oo) and does not depend on the "hierarchy" 
system for X operators. In the limit of infinite spatial 
dimensions, these analytical schemes allow u^to build 
a self-consistent Kadanoff-Baym type theoryO for the 
Hubbard model and some analytical results are given 
for simple approximations. The Falicov-Kimball model 
is also considered as an exactly soluble limit of Hubbard 
model. 



II. PERTURBATION THEORY IN TERMS OF 
ELECTRON HOPPING 

We consider the lattice electronic system that can be 
described by the following statistical operator: 



(2.1) 



a{(3)=TexpLfdT fdr' J^t^r ~ r')aUr)a,Ar') \ 
[ ^^'^ J 



where 



^0 — ^ Hi 



(2.2) 



is a sum of the single-site contributions and for the Hub- 
bard model we must put 

tf^{T~T')^U,S{T-T'). (2.3) 



In addition for the Falicov-Kimball model we must put 



It is supposed that we know eigenvalues and eigen- 
states of the zero-order Hamiltonian (2.2), 



H,\z,p) = Xp\z,p) (2.5) 
and one can introduce Hubbard operators 

Xp^\^,p){^,q\ (2.6) 
in terms of which zero-order Hamiltonian is diagonal 



(2.7) 



For the Hubbard model we have four states \i,p) = 
\i,n,T,nii): \i,0) = |i,0,0) (empty site), |i,2) = |i,l,l) 
(double occupied site), \i, t) = |«, 1, 0) and \i, |) = \i, 0, 1) 
(sites with spin-up and spin-down electrons) with ener- 
gies 

Ao = 0, \2 = U -2n, Xi = h- fi, A| = -ft. - fi. (2.8) 

The connection between the electron operators and the 
Hubbard operators is the following: 



X 
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aX, 



S2 



(2.9) 



Our aim is to calculate the grand canonical potential 
functional 



r! = -ilnSpp = r!o-^ln(a(/3))o, 
f^o = -^lnSpe^^^«, 



(2.10) 



single-electron Green functions 



G^Ar - r') = {Tal{T)a,^{T')) = ^'^^ (2.11) 



and mean values 



1 dO. 

i 

n — n-^ -\- riy] m — — n^, 



(2.12) 



where fXr^ = jj, + ah is a chemical potential for the 
electrons with spin a. Here, (...) = (1/Z) Sp(. . . p), 
Z = Sp (0, or in interacting representation 
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(^(/9)>o 



.am, 



.am 



Oc, 



(2.13) 



where (. . .)o = (l/Zo) Sp(. . . e"^^"); = Sp e 

We expand the scattering matrix a{P) in Eq. {2A) into 
the series in terms of electron hopping and for {a{P))o we 
obtain a series of terms that are products of the hopping 
integrals and averages of the electron creation and anni- 
hilation operators or, using Eq. (2J), Hubbard operators 
that will be calculated with the use of the corresponding 
Wick's theorem. 

Wick's theorem for Hubbard operators was formulated 
in Rcf. H (see also Ref. ^ and references therein) . For the 
Hubbard model we can define four diagonal Hubbard op- 
erators X'P'P [p = 0, 2, J,, t) which are of bosonic type, four 
annihilation X^^ , X^^ , X'^'^ , X^'^ and four conjugated 
creation fermionic operators, and two annihilation X'^\ 
and two conjugated creation bosonic operators. The 
algebra of X operators is defined by the multiplication 
rule 



xrx. 



PI _ 



^snX 



rq 



(2.14) 
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the conserving condition 



and the commutation relations 



(2.15) 



(2.16) 



where one must use anticommutator when both oper- 
ators are of the fermionic type and commutator in all 
other cases. So, commutator or anticommutator of two 
Hubbard operators is not a c number but a new Hub- 
bard operator. Then the average of a T products of X 
operators can be evaluated by the consecutive pairing, 
while taking into account standard permutation rules for 
bosonic and fermionic operators, of all off-diagonal Hub- 
bard operators X^'^ according to the rule (Wick's theo- 
rem) 



xr{T,)xi'^ (r) = -<5o.5p,(T-n)[xr(Ti),xr(n)]± 

(2.17) 

until we get the product of the diagonal Hubbard op- 
erators only. Here we introduce the zero-order Green's 
function 



9pq{ 



■ri) 



(2.18) 



±n±(Ap,) - 1 T < Ti ' 



where \pq = \p — \q and n±(A) = jftjjj a-nd its Fourier 
transform is equal 



(2.19) 



In particular, for the Hubbard model one can introduce 
the following pairings: 



+g,,{T-n){xPin) + xr{ri))}, 



a^^{Tl)aUT) = -%/,(r- ri)Xr(Ti), 



aUri)aUr) = <5.,/.(r - n) • a • (n), (2.20) 



a'hiTi)Xj'^{T) = Sijgasir - Ti)aJ^(Ti), 



I " — I 

ai„{Ti)XJ''{T) = -Sijgcrg{r - Ti)aia{Ti), 



a,„{Ti)Xf{T) = -Stjg2o{T - Ti) • • a^(Ti), 



where 



fallen) = gaoiuJn) " 52(7 (w„) 
= -Ugao{ljJn)g2s{'^n)- 



(2.21) 



Applying such pairing procedure to the expansion of 
(<5'(/3))o we get the following diagrammatic representa- 
tion: 



exp-! 




(2.22) 



where arrows denote the zero-order Green's functions ( ^.19| ), wavy lines denote hopping integrals and □, . . . stay for 
some complicated "n vertices" , which for such type perturbation expansion are an irreducible many-particle single- 
site Green's functions calculated with single-site Hamiltonian (2.7). Each vertex (Green's function) is multiplied by a 
diagonal Hubbard operator denoted by a circle and one gets an expression with averages of the products of diagonal 
Hubbard operators. 

For the Falicov-Kimball model expression ( 2.22| ) reduces and contains only single loop contributions 



(^(/3)>o 



where 



exp 




(2.23) 



p- 
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— Y and by introducing pseudospin variables 



= i {Pj^ — P~ ) one can transform the Falic ov-K imball model into an Ising-type model with th e eff ective multisite 
retarded pseudospin interactions. Expression (2.23) can be obtained from the statistical operator (2.1) by performing 
partial averaging over fermionic variables, which gives an effective statistical operator for pseudospins (ions). 

So, after applying Wick's theorem our problem splits into two problems: (i) calculation of the irreducible many- 
particle Green's functions (vertices) in order to construct expression (2.22) and (ii) calculation of the averages of the 
products of diagonal Hubbard operators and summing up the resulting series. 
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III. IRREDUCIBLE MANY-PARTICLE GREEN'S FUNCTIONS 



For the Hubbard model by applying the Wick's theorem for X operators one gets for two-vertex 

® = 5.oK)(Xr + + 92-A^n){Xr + Xf), 

for four-vertex 



ui^ia _ A (4) / \ 

-'^j°°5(To(t^n).9o-o(^^n+m) {U + C^^520 (^^ri+ri'+m)) 5ff0 (^^n' )35-0 (^^n'+m) 

n-\-n' -i-m 

+-^i^'^.925-(w„)52ff(^n+m) {U - U'^ gas{^n-7i')) 99o{^n')gso{^n' +m), 



(3.1) 



(3.2) 



and so on. Expressions ( |3.l| ) and (3.2) and for the ver- 
tices of higher order possess one significant feature. They 
decompose into four terms with different diagonal Hub- 
bard operators X^'p, which project our single-site prob- 
lem on certain "vacuum" states (subspaces), and zero- 
order Green's functions, which describe all possible exci- 
tations and scattering processes around these "vacuum" 
states: i.e., creation and annihilation of single electrons 
and of the doublon (pair of electrons with opposite spins) 
for subspaces p = Q and p = 2 and creation and annihila- 
tion of single electrons with appropriate spin orientation 
and of the magnon (spin flip) for subspaces p =| and 
P=i- 



In compact form expressions (3.1 
written as 



and (3.2) can be 



p 



Xfgaip) (Wn) 



and 



(3.3) 



(3.4) 



X C^(T(t(p) {^n,^l \^m)gs(p) {^l)9s(p) {^l+m), 



where 



ga(p){^n) — 

Here 



gao{LOn) for p= 0,(7 
325- (^n) forp=:CT,2 



(3.5) 



U ± 1/^320 (t^n+/+m) for p = 0, 2 
U ± U^gaai^^n^l) for p = (T, 



Ucra(p){^n,^^l\^^m) = C^ctct(p) (w; , CJ„ |cJm) 



(3.6) 



is a renormalized Coulombic interaction in t he su bsp aces. 
In diagrammatic notations expressions (3^) or (3^) can 
be represented as 




for p — (T,a 



(3.7) 



where dots denote Coulombic correlation energy U — 
•^2 + '^0 ^ -^T ^ '^i dashed arrows denote bosonic zero- 
order Green's functions: doublon (720 (wm) or magnon 

ffcrS-(Wm)- 

For six- vertex one can get 



X X! ^f^5cr(p) (w„)5cr(p) (^m )9s{p) (W„2 )gg(p) {i^ns )9a(p) (w„4 )ga(p) (^n^ ) 
P 
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— Ua9(p) {^n , I ljJ„2 )9s(p) {^n+n2 -ns ) ^^crff (p) ('^"l i "^"4 l^ns -"4 ) 

— Ucra(p) {^n, |aj„^„„3)55(p) (ti'n+n4-ri3)t^<Tff(p) (^ni j '^n2 l^ns-ng) 

+ '^(Tffff (p) ("^n : '^Tii , j ^^ria , , ) | , 

where 

^ a3S(p) {^n , IjJni , , , , OJn^ ) 

^ f ±f^^ (520(^n+ri2) - 320(w„+„4)) (^20 (t^m +ri3 ) ^ g2o(w„i+„5)) for p 

— 714 )) for p 



(3.9) 



0,2 
a, a 



In expression ( |3.8| ) the contributions of the first four 
terms in braces can be presented by the following dia- 
grams: 



(3.10) 



with the internal vertices of the same type as in Eq. ( p.7D , 
whereas the contribution of the last term can be pre- 
sented diagrammatically as 



>V-< 

So, we can introduce primitive vertices 

X>- X 



(3.11) 



(3.12) 



by wh ich one can construct all n vertices in expansion 
( ^.22 ) according to the following rules: 



1. n vertices are constructed by the diagonal Hubbard 
operator X^p and zero-order fermioni c and bosonic 
lines connected by primitive vertices ( 3.12|) specific 
for each subspace p. 

2. External lines of n vertices must be of the fermionic 
type. 

3. Diagrams with the loops formed by zero-order 
fermionic and bosonic Green's functions are not al- 
lowed because they are already included into the 
formalism, e.g., gives r'^ 

For n vertices of higher order a new primitive vertices 
can appear but we do not check this due to the rapid in- 
crease of t he al ge braic calcul ation s with the increase of n. 
Diagrams (3/7), ( 3.1C ), and ( 3.11 ) topologically are trun- 
cated Bethe-lattices constructed by the primitive vertices 
( ^l2| ) and can be treated as some generalization of the 



Hubbard starsLlil-j in the thermodynamical perturbation 
theory. 

It should be noted that each n vertex contains Coulom- 
bic interaction U as in primitive vertices ( 3.12| ) (denoted 
by dots) as in th e denominators of the zero-order Green's 
functions (2.19). In the C/ — s- oo limit, each term in the 
expressions for n vertices can diverge but total vertex 
possesses finite U oo limit when diagrammatic series 
of Rcf. 1^ are reproduced. 

The second problem of calculation of the averages of 
diagonal X operators is more complicated. One of the 
ways to solve it is to use semi-invariant (cumulant) ex- 
pansions as was done in Refs. ^ and |^ for the U = oo 
limit. Another way is to consider the d — oo limit where 
new simplifications appear. 



IV. DYNAMICAL MEAN-FIELD THEORY 

Within the frames of the considered perturbation the- 
ory in ter ms of electron hopping a single-electron Green's 
function (2.11) can be presented in a form 



(4.1) 



where we introduce an irreducible part [ujn , k) of 
Green's function which, in general, is not local. In the 
case of infinite dimensions d oo one should scale the 
hopping integral according to 



t. 



Vd 



(4.2) 



in order to obtain finite density-of-jS±ates and it was 
shown by Metzner in his pioneer worked that in this limit 
the irreducible part become local 

Sjjo-(T - r') = (5ySo-(T - t') or E„{ujn,k) ^ E^{uJn) 

(4.3) 
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and such a site-diacppnal function, as it was shown by 
Brandt and Mielschjj can be calculat ed b y mapping the 
infinite-dimensional lattice problem (2.1) with tfj{T — 
t') = -^tij5{T — t') on the atomic model with auxiliary 
the Kadanoff-Baym field 



t?^(T-r') = %J.(r-r'), 



(4.4) 



which has to be self-consistently determined from the 
condition that the same function So.(wn) defines Green's 
functions for the lattice and atomic limit. The self- 
consistent set of equations for 'E.„{ujn) and Jo-(w„) (e.g., 
see Ref. ^ and references therein) is the following: 

1-^1 _ 1 

= Gi'^)(c.„,{J,(c^„)}), (4.5) 

where ci"^ (w„, { Jo-(cj„)}) is a Green's function for 
atomic limit (4.4). 

Grand canonical potential for lattice is connected with 
the one for atomic limit by the expressioncl 



lnGi'^)(c.„)~^Eli^^-(^"'^)f 

k J 



(4.6) 



On the other hand, we can write for the grand canon- 
ical potential for atomic limit f2a the same expansion as 
in Eq. ( 2.22| ) but now we have averages of the products 
of diagonal X operators at the same site. According to 
Eq. (2.14) we can multiply them and reduce their product 
to a single X operator th at can be taken outside of the 
brackets and exponent in (2.22) and its average is equal 
to (XPP)o = J: ^-^x . Finally, for the grand canonical 
potential for atomic limit we get 



1 

7" 



(4.7) 



where 




+ 




(4.8) 



DO 



r 



are the "grand canonical potentials" for the subspaces. 

Now we can find single-electron Green's function for 
atomic limit by 



sn 



T — T 



SJa{T 



where 



Go-(p)(t- r') 



sn 



(4.9) 



(4.10) 



are single-electron Green functions for the subspaces 
characterized by the "statistical weights" 



g-0s:2(p) 



(4.11) 



and our single-site atomic problem exactly (naturally) 
splits into four subspaces p = 0, 2, J,, f. 



We can introduce irreducible parts of Green's functions 
in subspaces So.(p)(w„) by 



Gct(p)(Wji) — 



1 



where 



(4.12) 



(4.13) 



According to the rules of the introduced diagrammatic 
technique, n vertices are ter mina ted by the fermionic 
Green's functions [see (It]), (Fio|) , and (|3ll| )] and this 
allows us to write a Dyson equation for the irreducible 
parts and to introduce a self-energy in subspaces 



ct(p)^ 



-cr(p) 



(4.14) 



where self-energy I]cr(p)(^n) depends on the hopping in- 
tegral Ja' (wn' ) only through quantities 



*CT'(p)(^n') = GCT'(p)(a;„') — ^o.'(p)(^n') 



(4.15) 



"(t'(p) 



{Un' )Ja'{uJn') + ^a'{p){t^n')Ja' (w„' ) H } 
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It should be noted, that the total self-energy of the 
atomic problem is connected with the total irreducible 
part by the expression 



(4.16) 



and it has no direct connection with the self-energies in 
the subspaces. 



The fermionic zero-order Green's function (B.5) can be 
also represented in the following form 



1 



9'^(p) - . (0) 



where 



ct(p) 



dXr, 



a(p) 



for p = 0, (T 



dficr \ 1 for p = 2, (7 



(4.17) 



(4.18) 



is an occupation of the st ate \p ) by the electron with spin 
a, and Green's function ( [1.12| ) can be written as 

Ga[p){l-On) — T^J) • 

(4.19) 

Now, one can reconstruct expressions for the grand 
canonical potentials ^{p) in subspaces from the known 
structure of Green's functions. To do this, we scale hop- 
ping integral 

J^(a;„) ^ aJA^n), a £ [0, 1], (4.20) 

which allows to define the grand canonical potential as 
1 

^(p) = \ + J da^'^J„{LJn)G^(^p){u;n,a) (4.21) 

and after some transformations one can get 
1 



■El- 



(4.22) 



where 



(p) 



1 

^E Jda^aip) {cJn,a) 



da 



(4.23) 



is some functional, such that its functional derivative 
with respect to 5* produces self-energy: 



S<i> 



(p) 



'^^<T(p)('^n) 



rip)\ 



(4.24) 



So, if one can find or construct self-energy Scr(p) i'-^n) he 
can find Green's functions and grand canonical potentials 
for subspaces and, according to Eqs. (17) and ( [4. 9] ), solve 
atomic problems. 



Sta rting from the grand canonica l pot ential (4.7) and 
( 4.22 ) one can get for mean values ( 2.12| ), 



J^WpU^ip), 



(4.25) 



"-(P) = ^aip) ^^2^ 



(0) 



[Gcr(p){i^n) - ^a{p){^n)\ 



9$ 



(p) 



where in the last term the part ial derivative is taken over 
the Ha not in the chai ns (|4.15| ). The second term in the 
right-hand side of Eq. ( 4.25| ) can be represented diagram- 
matically as 



(4.26) 

and the first contributions into the last term are following 




(4.27) 

where double lines denote quantities a(p){^n)- Loop 
is connected with the superconducting or magnon 
susceptibilities for subspaces p = 0, 2 or p = cr, (t, respec- 
tively. 

For the single atom [Jo-((jj„) = 0] we have $(p) = 0, 

Ga(p){^n) = '^a(p){^n) = 5cr(p)(^n); a-^d 
p n p 

but in the general case [Jo-(a;„) ^ 0] we cannot prove that 
the sum rule 



(4.29) 



is fulfilled. 



A. Falicov-Kimball model 

For the Falicov-Kim ball model J|(cj„) = and accord- 



ing to Eqs. (BJ) and (3 



5^T(p)('^") = 0; ■=T(p)('^") = ffT(p)('^") (4-30) 



and 



^(p) = \ - 4 E 1^ ~ -^T ('^«)5t(p) M] , (4.31) 



(4.32) 
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(4.33) 



which immediately gives results of Ref. || (see also 
Ref. H). 

For the Hubbard model there are no ex act expr ession 
f or se lf-energy but the set of Eqs. ( 4.12| ), ( 4.14 ), and 
(4.22) allows one to construct different self-consistent ap- 
proximations. 



B. Alloy-analogy approximation 

The simplest approximation, which can be done, is to 
put 



which gives 



and 



S<T(p)(Wn) = 



Scr(p)(Wn) = ga{p)i^n) 



(4.34) 
(4.35) 



^(p) = -^P " 4 X! ^^'^ ^ J<7{i^n)9cr{p){^n)] (4.36) 

and for the Green's function for the atomic problem one 
can obtain a two-pole expression 



Wo + Wa 



W2 + Wg 



iuJn — AcrO — Ja{uJn) i^n — ^2S — J<j{^n) 

(4.37) 



of the alloy-analogy solution for the Hubbard model, 
which is a zero-order approximation within the consid- 
ered approach and is exact for the Falic ov-K imball model. 
For this approximation, mean values ( ^.12| ) are equal to 



Wq + Wa W2+ Wg 

p^-<y v-'w g/3A.o +1 ^ eP^-^^ +1 



^'y = \Yl G'a'' (Wn) +W2+Wa 
n 

^i^G(-)(^„) 



(4.38) 



and, for some values of the chemical potential, they can 
get unphysical values: negative or greater then one. 



C. Hartree-Fock approximation 

The next possible approximation is to take into ac- 
count the contribution from diagram (4.26) and to con- 
struct the equation for the self-energy in the following 
form: 



M^'^Y^ U^sf^p) [iUn' ) , (4.39) 

n' 

which, together with the expression for mean values 



^(t(p) 



V(p) 



(4.40) 



- X + ^ tanh ^ [U n^^p) - 



<y[p) 2 ' 2 2 

n' 

gives for the Green's function in the subspaces expression 
in the Hartree-Fock approximation: 



1 



(4.41) 



Now, grand canonical potentials in the subspaces are 
equal 

^(p) = K - IsY^'^i^ - Ja{l^n)^a{p) (w„)] 
ncr 

-U (n.(p) - - (4.42) 



and for the Green's function for the atomic problem (4.9) 
one can obtain a four-pole structure 



G('^)(c.„)-E: 



(4.43) 



Expr ession ( 4.43| ), in contrast to the alloy- analogy solu- 
tion (4.37), possesses the correct Hartree-Fock limit for 
small Coulombic interaction J7 <C t: 



1 



(4.44) 



when Wp « I and n^(^p) n„ = ^ Y.n G f iuj„). On the 
other hand, in the same way as an alloy-analogy solu- 
tion, it describes the metal-insulator transition with the 
change of U. 

In Fig. 0the frequency distribution of the total spectral 
weight function 



(4.45) 



as well as cont ribut ions into it from the subspaces [sepa- 
rate terms in ( fl:.43|) ] are presented for the different elec- 
tron concentration (chemical potential) values. One can 
see, that the spectral weight function contains two peaks, 
which correspond to the two Hubbard bands. Each band 
is formed by the two close peaks: p = and a for the 
lower Hubbard band and p — 2 and a for the upper one, 
with weights Wp Eq. (4.11). The main contributions come 
(see Fig. g) from the subspaces p = for the low electron 
concentrations {n<^,fj,<0),p — 2 for the low hole 
concentrations (2 — n < |, /i > [/) and p = cr,a for the 
intermediate values. For the small electron or hole con- 
centrations, the Green's function for the atomic problem 
(4.43) possesses correct Hartree-Fock limits too. 
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FIG. 



1. Spectral weight function pa{ijj) ( 4.45 ): 
n = 1; (b) ^ = -1, n = 0.07; (c) = 0.01, 



total and for each subspace. 
n = 0.72; (d) /i = -0.01, n 



for the different chemical potential values: (a) 
= 0.66 (f/ = 4, r = 0.2). 




FIG. 2. Statistical weights of the subspaces Wp ( 4.11 ) as 
functions of the electron concentration (U = 4, T — 0.2). 

Such four-pole structure of the single-electron Green's 
function can be obtained also for the one-dimensional 
chain with the N = 2 periodic boundary condition 
(see the Appendix), which is equivalent to-the two-site 
problem considered by Harris and LangcEj Here, two 
poles correspond to the noninteracting electrons or holes, 
which hope over the empty sites, and give the main con- 



tribution for small concentrations. The other two poles 
give the main contribution close to half-filling and corre- 
spond to the hopping of the strongly-correlated electrons 
over the resonating valence bond (RVB) states. 

So, one can suppose that the Hubbard model describes 
strongly-correlated electronic systems that contain four 
components (subspaces). Subspaces p = and p — 2 
describe the Fermi-liquid component (electron and hole, 
respectively) which is dominant for the small electron and 
hole concentrations, when the chemical potential is close 
to the bottom of the lower band and top of the upper one. 
On the other hand, subspaces p =1 and | describe the 
non- Fermi- liquid (strongly correlated, e.g., RVB) compo- 
nent, which is dominant close to half-filling. Within the 
considered Hartree-Fock approximation, at n « | and 
2 — n w |, we have transition between these two regimes: 
Fermi liquid and non-Fermi liquid. It reminds us the 
known properties of the high-Tj, compounds, where for 
the nondoped case (n = 1) compounds are in the an- 
tiferroelectric dielectric state, then for small doping the 
non-Fermi-liquid behavior is observed (underdoped case 
n < 1) and after some optimal doping value, the proper- 
ties of the compound sharply change from the non- Fermi 
to the Fermi liquid (overdoped case). 

The results presented in Figs. |^ and || are obtained 
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for relatively high temperature. With the temperature 
decrease, on the one hand, the transition between the 
Fermi and non-Fermi liquid becomes sharp and, on the 
other hand, for some che mical potential values there can 
be three solutions of Eq. (4.4C) with two of them corre- 
sponding to the phase-separated states. The considera- 
tion of the phase separation in the Hubbard model is not 
a topic of this paper and will be the subject of further 
investigations. 



D. Beyond the Hartree-Fock approximation 

Se lf-ene rgy in the Hartree-Fock approximation [see 
Eq. (4.41)] describes some self-consistent shift of the ini- 
tial energy levels and does not depend on the frequency. 
All other improvements of the expression for self-energy 
add the frequency dependent contributions. To see this, 
let us consider the contribution into the mean values from 



U 




the first diagram in Eq. ( 4.27 ). This diagram originates 
from the following skeletal diagram 



(4.46) 



in the diagrammatic expansion for functional *&(p). On 
the other hand, such a skeletal diagram produces addi- 
tional contribution into the self-energy 



which is frequency dependent. Also, in order to get a 
self-consistent set of equations, we introduce renormal- 
ized bosonic Green's functions 



£'20 (^ 



1 



1 



ioJm — A20 ioJm — Actct 

A20 = A20 + c/^Xl*'^(p)('^")' 



Finally, for the Green's function ( 4.19 ) we get the gen- 
eral representation 



1 



(4.49) 



where the Hartree-Fock contribution U ng(^p^ is extracted 

and SCT(p)(a;„) is a frequency dependent part of the self- 
energy, which within the considered approximation is 
equal 



(4.50) 



where 



Sa{p)i^7i) = ±— ^D„<j(p)(a;„,cj„/)^'ff(p)(cj„/) (4.51) 



D2o{'^n+n') for p = 0,2 
D„9 (ujn-n- ) iov p = a,a 



(4.52) 



and 

Now, mean values ( [4.25| ) are equal 

n 

nn' 

and for the grand canonical potentials in the subspaces 
we obtain the following expressions: 

^(p) = -^p ^ 4 X!^'^ ^ J^(w„)S^(p)(w„)] (4.54) 
C^-D|o(w„+„05^o^(^^«+«')] 



X*a(p) (W„)5'ff(p) [uJn' 



for p — 0,2, and 



^(p) ^ \ ~ ^^^^'^{^^ (^«)S^(p) i^^n)) (4.55) 



/3' 



^ ^ (7 (1 ± C/i?2^(^„_„0ff^i(^n-n')) 



X*c7(p) ('^n)*ff(p) (^ri 



for p = a,(T. 

In order to analyze the structure of the poles in 
Eq. ( 4.49 ), an analytical continuation of the expression 
for Eo.(-p)(a;„) from the imaginary axis to the real one 
should be done. To do it, we use the well-known identity 



-T 



(3 ^ iujn — A 



- ±n±(A), 



(4.56) 



which follows from Eq. ( 2.18 ), and analytical properties 
of the Green's function 



TT Z — Ll! 



(4.57) 



Green's functions in the subspaces G„^p){z), irreducible 
parts 2o.(p)(-2^); and dynamical mean-field Ja{z) all pos- 
sess the same analytical properties. Finally, for S'o-(p)(z) 
we get the following expressions: 



Q A-^-uf^ 5^^(p)(^>-jO+) 

TT 7 Z + UJ-X20 

— OQ 

±n_ (A20) C/*ff(p) (A2o~z) (4.58) 
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for subspaces p = 0,2 and 



V. SUMMARY 



Sa(p){z)^±—V dujn+{uj) 



Z — U) — Actct 

n^{\„s) + l ?7*ff(p)(z- A^ff) (4.59) 



for p — g, a . An alytic al continuation of expressions 
(4.53), (4.54), and (4.55) can be done in the same way. 
One can see, that contributions (4.58) and ( 1.5S| ) diverge 
when A20 — and \as — 0, respectively, which is an 
unphysical result. 

So, we cannot include into the consideration only one 
contribution from diagram (4.47) but must sum up all 



diagrams of the following type 



(4.60) 



which gives an expression free from the above-mentioned 
divergences 



K)-;:^(i + ^;'(p)K) 



u 



l + 5'^,(p)(Wn 



1 + 'S'cr(p)('^n) 

(4.61) 



where S„(^p){ijJn) is defined above and 



/3 



(4.62) 



Sa(p){^n) = j^D^g(^p-){uJn,UJn')'^a{p){i^n') 
n' 

I 1 + S's(p)(w„') 



Such diagram resummation m ust b e also done in the ex- 
pression for the mean values (4.53), where the last term 
must be replaced by 



±- 



1 



^ ^'^'-(P) ('^" ' (^n)*ff (p) (w„' ) (4.63) 



^ + Kip)M , i + 5;w(t^n') 



1 + '5'<j(p)(w„) 1 + S'ct(p)(w„') 



- 1 



Besides diagram ( 4.47 ), there are a lot of other dia- 
grams that diverge and need additional resummation of 
the diagrammatic series. But now it is difficult to clear 
out what types of diagrams are leading in different case, 
which calls for additional investigation. But it is, ob- 
viously, that such contributions will shift the boundary 
between the Fermi and non-Fermi-liquid behavior. 



A finite-temperature perturbation theory scheme in 
terms of electron hopping, which is based on the Wick 
theorem for Hubbard operators and is valid for arbitrary 
values oi U {U < 00) has been developed for Hubbard- 
type models. Diagrammatic series contain single-site ver- 
tices, which are irreducible many-particle Green's func- 
tions for unperturbated single-site Hamiltonian, con- 
nected by hopping lines. Applying the Wick theorem 
for Hubbard operators has allowed us to calculate these 
vertices and it is shown that for each vertex the problem 
splits into subspaces with "vacuum states" determined by 
the diagonal (projection) operators and only excitations 
around these "vacuum states" are allowed. The vertices 
possess a finite U 00 limit wkep diagrammatic series 
of the strong-coupling approachHEl are reproduced. The 
rules to construct diagrams by the primitive vertices are 
proposed. 

In the limit of infinite spatial dimensions the total aux- 
iliary single-site problem exactly (naturally) splits into 
subspaces (four for Hubbard model) and a considered an- 
alytical scheme allows to build a self-consistent Kadanoff- 
Baym-type theory for the Hubbard model. Some ana- 
lytical results are given for simple approximations: an 
alloy-analogy approximation, when two-pole structure 
for Green's function is obtained, which is exact for the 
Falicov-Kimball model, and the Hartrec-Fock-type ap- 
proximation, which results in the four-pole structure for 
the Green's function. Expanding beyond the Hartree- 
Fock approximation calls for considering of frequency de- 
pendent contributions into the self-energy and resumma- 
tion of the diagrammatic series. 



In general, the expression 



(5.1) 



+ - Ung(^p) - So.(p)(wn) - J(T(a;„) 



gives an exact four-pole structure for the single-electron 
Green's function of the effective atomic problem. In 



Eq. (4.14) zero-order Green's functions (3.5) are the same 
for the subspaces p = Q,a and p = 2, ct, respectively, and 
correspond to the two-pole solution of the one-site prob- 
lem without hopping. Switching on of the electron hop- 
ping splits these two poles and the value of splitting is 
determined by the values of the self-energy parts in the 
subspaces, which describe the contributions from the dif- 
ferent scattering processes. Alloy-analogy approximation 
neglects by the such scattering processes iYi„(^p-^{uJn) = 0) 
which res ults in the two-pole structure for the Green's 
functions (4.37). But, in general. Green's functions pos- 
sess four-pole stru cture and even the Hartree-Fock ap- 
proximation (4.43) clearly shows it. 
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It should be noted that the four-pole structure of the 
Green's function for the atomic problem might not re- 
sult in the four bands of the spectral weight function 
(see Fig. 0). The presented consideration allows us to 
suppose that each pole describes contributions from the 
different components (subspaces) of the electronic sys- 
tem: Fermi liquid (subspaces p — 0,2) and non-Fermi 
liquid {p =t,i), and for small electron and hole concen- 
trations (n < I and 2 — n < |) the Fermi-liquid compo- 
nent gives the main contribution ( "overdoped regime" of 
high-Tc's), whereas in other cases the non- Fermi liquid 
one ( "underdoped regime" ) . 



ACKNOWLED GMENTS 



I am grateful to Professor I. V. Stasyuk for the very use- 
ful and stimulating discussions. Also I acknowledge the 
hospitality of the Institute of Physics of the M. Curie 
Sklodowska University of Lublin and thank Professor 
K.I. Wysokihski for discussions. 



APPENDIX: TWO-SITE PROBLEM 



Let us consider an infinite one-dimensional chain with 
the N = 2 periodic boundary condition. Mathemati- 
cally it is equivalent |-to the two-site problem considered 
by Harris and Langejl3 but now we can introduce the lat- 
tice Fourier transformation, with two wave-vector values 
in the first Brillouin zone q = and q — n, and perform 
all calculations for the grand canonical ensemble. The 
Hamiltonian of the model is the following: 



1=1.2 \ a I 



(Al) 



We can introduce the Fourier transform of the electron 
hopping 



< for g = 
tn ^ tcosq ^ < , r 
^ ^ ' —t tor q = n 



(A2) 



and our aim is to calculate the single-electron Green's 
function 

, ^ _ ( Giia{uJn) + Gi2aii^n) for q = . 

^-^^^h) - <^ GnA^n) - Gi2aM for 9 = ^ ' ^^^^ 

where G,,,(t - r') = (TaUT)a,AT') 

The initial basis of states contains 16 many-electron 
two-site states \pi,P2), where pi — {n,i|n^}, i.e., 



Il) = 


0,0), 








|2) = 


1,0) 


= a\,\l). 






|3> = 


o,i> 


= 4jl). 






|4) = 


T,o) 


= al||l), 






|5>- 


o,T> 


= 4t|1)' 






|6) = 


i,i> 


= 4J2) = 


-alj3). 




|7) = 


T,T) 


= 4t|4) = 


-al^\5), 




|8>- 


2,0) 


= «1t|2> = 


-«u|4), 




|9) = 


T,i> 


= aJtl^) = 


-4J4), 




|io>- 


i,T> 


= 4t|2> = 


-alj5). 




|11) = 


0,2) 


= 4t|3> = 


-4j5), 




|12) = 


2,i) 


altl^) 


alj9)=- 


4j8), 


|13) = 


i,2) 


= al^\^) = 


atjU) = 


-4J10), 


|14)- 


2,T) 


-«li|7> = 


4tI8> = - 


4^110), 


|15) = 


T,2) 


= 4il7) = 


4tI9> = - 


-4^111), 


|16) ^ 


2,2) 


-4t|13) = 


= -4t|12) 


= 4 J 15) 



(A4) 



and one can introduce Hubbard operators X^'^ — \p){q\ 
acting in the space of these states. Now, electron creation 
operators can be presented in the following form: 



X 



4,1 



X 



7,5 



^8,2^^9,3 
_|_^12,6 _ j^l4,10 _ j^l5,ll 
2,1 _ -^6,3 _ j^8,4 _ j^l0,5 



(A5) 



^5,1^^7,4^^10,2^^^11,3 



16,12 



= X 



3,1 



x" 



-X^^-'^ -X 



~X" 

13,10 



4 _ ^11,5 
^15,7 



X' 



By transformations 



|2) |4) |12) |14) 
|3) |5) |13) |15) 



(A6) 



V2 V2 



1 
V2 



|2) |4) |12) |14) 
|3) |5) |13) |15) 




sin (/) ^ 
COS0 ^ 
sin (A 4= cos 6-^0 



1 

± 

t 
t 

%/2 



sm ( 



C0S( 



t 

%/2 



sin( 



V2 





1 

V2 



|9) 

V|ii>/ 



where sin2(/i(t) = 
be diagonalized 



2t 



^[72/4+4*2 



(AT) 

, the Hamiltonian ( |Al[ ) can 
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H 



(A8) 



Here 



Ai 



A2 — Aj 

h = h 



"10 
As 



Ar, = 



^T2 = = 



0, 

-H + t, 

U + J- 2/1, 
-J-2/i, 
[/ - 2/i, 

[/ - + <, 

2t/-4/i 



(A9) 



are eigenvalues and 



4*2 



^C/2/4 + 4t2 + [//2 U 



(t/>i)- (AlO) 



Finally, with the use of the Wick theorem ( ^.17 ) for the 
Hubbard operators acting in the space o f eig enstates \p) , 
for the single-electron Green's function (|A3|) we obtain 



U! — t„ UJ + J + tn 



(All) 



+- 



UJ ■ 



U -tn OJ -U - J + tn' 



where 
4 



1 
Z 



1 

Z 



^-|3{2U-i^i) ^ g-/3(C/-3M+t) 



(A12) 



g/3(p-t) _,_ g/3(2p+J) g-/3(;7-2p+J) _^ g-/3(C/-3M-t) 



and Z = EpG"^^''. 

One can see, that Green's function (All) possesses 
a four-pole structure and the spectrum contains four 
"bands" grouped near the initial energy levels of the one- 
site problem and U . The distance between the centers 



of gravity of the grouped bands is equal to J [Eq. (AlC)] 
and is of the order of magnitude of the effective exchange 



interaction. It is obvious that the weights of the bands 
satisfy the sum rule 



Aim) + Mh) + Bi{tq) + B2m) ^ 1. 

The spectral weight function is equal 
p,{E) = ^Y.^G,{E-iQ+,t,) 



(A13) 



-[Ai{t)5{E~t)+Ai{-t)5{E + t) 



(A14) 

+A2{t) 5{E-U-t)+ A2{-t) S{E-U + 1) 
+Bi{t) d{E + J + t)+ Bi{~t) 6{E + J- t) 
+B2{t) S{E - U - J + t) 
+ B2{-t)S{E -U - J -t)], 

contains the same eight energies obtained by Harris and 
Langellj but with different weights and originates from 
the four poles (bands) of the Green's function (All) for 
the two-site problem. 

The nature of these peaks is clear from the ground- 
state properties of the model. At zero temperature, de- 
pending on the value of the chemical potential or elec- 
tron concentration, the ground states are the following 
{U :S> t): empty state n — < — <): 



|1> = |0,0), 
Ai =0, 



(A15) 



one-electron states n = ^ (— t < fi < t — J): 





V«2i 








^^2T 


"it 




A5 = 




-t, 



G^ico,q = 0) = 
Ga{uJ,q = tt) 



3/4 1/4 



(A16) 



UJ — t UJ ~ U — t 
1/2 i(l-Hsin20) i(l-sin20) 



UJ + t UJ + J — t UJ — U — J — 

two-electron state n — \ {t — J < fi < U + J — t): 

|9> = ^ cos (a\^al^ - "u^t) I^' 0) 
-^sin(/.(^a|;^ai^-a^^4^) |0,0), 



G„{ui,q) 



2/i- J, 

i(l-hsin20q) i(l -sin2(/),) 



UJ + J + t„ 



UJ -U - J + t„ 



(A17) 



three-electron states n = ^ {U + J — t < fi < U + t): 
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|13) = -^(aiT+a2T)|2,2), 

|15)--^(au+a2i)|2,2), 
Ajg = Ap^ = U — 'Sfi ~ t, 

Ga(LU,q = 0) - 



LU + J + t 

i(l + sm20) 1/2 



(A18) 



Ga{uj,q = tt) = 



oj-U-J + t oj-U-t' 
1/4 3/4 



tJ + t UJ - U + t' 
and four-electron state n — 2 {fj, > U + t): 

|T6) = |2,2), 
Ajg = 2C/ - 4^^, 

1 



uj ~ U ~ t„ 



For small electron (n « 0) or hole 



(A19) 



concen- 



trations we get Green's functions (A15) and (Al£), re- 
spectively, which describe hopping of the noninteracting 
particles over empty states. 

On the other hand, for the half-filled (symmetric) case 
n w 1, the ground state 19) is mainly a RVB-type state. 
Now, Green's function ( A17 ) possesses two-poles shifted 
by the value of the effective exchange interaction J from 
the one-site levels and describes the electron transfer over 
the RVB states. The weight of each pole depends on the 
hopping value, but its total contribution into the spectral 



weight function (A14) is equal to ^ as it should be for 



the symmetric case. 

For other cases the number and weights of the poles 
in the spectral weight function (A14) strongly depend on 
the electron concentration (chemical potential) and wave- 
vector values and contain contributions from the nonin- 
teracting electrons (holes) and the strongly hybridized 
RVB states. 
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